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ABSTRACT 

Broad Fe II emission is a prominent feature of the optical and ultraviolet spectra of quasars. We report 
on a systematical investigation of optical Fe II emission in a large sample of 4037 z < 0.8 quasars selected 
from the Sloan Digital Sky Survey. We have developed and tested a detailed line-fitting technique, taking into 
account the complex continuum and narrow and broad emission-line spectrum. Our primary goal is to quantify 
the velocity broadening and velocity shift of the Fe II spectrum in order to constrain the location of the Fe 
ll-emitting region and its relation to the broad-line region. We find that the majority of quasars show Fe II 
emission that is redshifted, typically by ~ 400 km s"' but up to 2000 km s~', with respect to the systemic 
velocity of the narrow-line region or of the conventional broad-line region as traced by the H/3 line. Moreover, 
the line width of Fe II is significantly narrower than that of the broad component of H/3. We show that the 
magnitude of the Fe II redshift correlates inversely with the Eddington ratio, and that there is a tendency for 
sources with redshifted Fe II emission to show red asymmetry in the H/3 line. These characteristics strongly 
suggest that Fe II originates from a location different from, and most likely exterior to, the region that produces 
most of H/3. The Fe ll-emitting zone traces a portion of the broad-line region of intermediate velocities whose 
dynamics may be dominated by infall. 

Subject headings: galaxies: nuclei — (galaxies:) quasars: emission lines — (galaxies:) quasars: general — 
galaxies: Seyfert — line: profiles 



1. EVITRODUCTION 

Fe emission contributes significantly to the optical and ul- 
traviolet (UV) spectra of most active galactic nuclei (AGNs), 
both in terms of wavelength coverage and flux. The proper- 
ties of the Fe ll-emitting clouds may provide important clues 
to the underlying physics in the broad-line region (BLR). 
First, Fe II emission can be used to constrain the cover- 
ing factor of BLR clouds from energy budget considera- 
tions. Second, the ratio of the equivalent width (EW) of Fe 
II to that of H/3bc strongly varies with statistical measures of 
AGN correlations, such as the so-called Eigenvector 1 derived 
from principal component analysis, which is believed to be 
driven by some fundamental p roper t y such as mass accretion 
rate (e.g., iB oroson & Green 1992t ISulentic et af] l2000allbL 
lMarziani"et aLi2001i , i2003ai) . Third, Fe abundance derived 
from Fe II emission can be used to study the cosmological 
evolution of AGNs and possib ly chemical enrichment of their 
hosts and environment (e.g.. iWHs et alj|1985l: IWheeler et al] 
fTosllDietrich et alJl200l l2003HMaiolino et al.ll2003h ." Care- 
ful measurement of the properties of Fe II emission in a large 
sample of AGNs will clearly have a considerable impact on 
our understanding of these systems. 

The origin of the optical/UV Fe emission has been hotly 
debated for more than two decades. Thousands of Fe emis- 
sion lines blend together to form a pseudo-continuum, which, 
when combined with Balmer continuum emission, results 
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in flie "small b lue bump" around 3000 A dGrandil [19821: 
IWills et aLlll985h . Previously, theoretical calculations of pho- 
toionized clouds in the BLR encountered difficulties repro- 
ducing the observed strength of strong Fe II emission, prompt- 
ing many authors to propose additional ph ysica l mechanisms 
(Netzer & Wifls 1983; Wifls et al. 198 51 iJoivl [19871 [19911: 
ICoUin-Souffrin et aljl988ciSigut & Pradhan'1998V Recently, 
a more sophisticated calculation by Baldwin et al. (2004) re- 
vealed that the predicted shape and EW of the 2200-2800A Fe 
II UV bump can only be made consistent with observed val- 
ues if either microturbulence of hundreds of km s"' or another 
collisionally excited component is included in the model. 

Despite its significance, current observations of Fe II emis- 
sion provide poor constraints on its origin. Some studies 
suggest that Fe II emission originates from the sam e region 
as the other broad emission line s. For example. [PhillipsI 
^I9ff), Boroso n & Gre^ ([19921) . iLaoretalJ (I1997H) . and 
[Wron-Cettv et al] (120041) obse rved similar line widths and 
profiles for Fe II and H/3bc, and lMaoz et alj (Il993h found that 
both Fe II emission and the Balmer continuum have compa- 
rable varia tion amplitudes. But there is a growing debate on 
this issue. Marziani et alj (l2003 b) found that H/3bc is system- 
atically broader than Fe II for sources with FWHM(H/3bc) > 
4000 km s"' . Recent studies of Fe II emission variability have 
shown that Fe II emission responds to variations in contin- 
uum flux but that th e variability amplitude of Fe II is not the 
same as that of H/3 ("Vestergaard & Peterson"2005' and refer- 
ences therein; Wang et al. 2005; Kuehn et al. 2008). In fact, 
the upper limit on the time lag between Fe II emission and the 
continuum exc eeds the lag of any other obs erved emission 
lines obtained dVestergaard & Peterson|[2005h . This implies 
that Fe II may be emitted from fu rther out in the BLR than 
any other broad emission line. Kuelmetal] ( l2008l) suggested 
that Fe II may be produced from a region betwe en the BLR 
and th e dust sublimation radius. The study of Matsuoka et alj 
(|2008|) . based on measurements of the O I and Ca II emis- 
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sion lines, supports the notion that Fe II emerges from the 
outer portion of the BLR. In a recent three-dimen sional spec- 
ti-oscopic study of Mrk 493. lPopovic et all (l2007l) found that 
the Fe II emission region is extensive and that the line width 
of Fe II is only 1/3 of that of H/3bc, leading them to suggest 
that Fe II emission originates in an intermediate-line region. 
If these findings can be confirmed, Fe II emission can be a 
probe of the intermediate-line region, which may be the tran- 
sition from the torus to the BLR and accretion disk. 

Quasar emission lines often exhibit considerable velocity 
shifts with respect to each other However, the relative ve- 
locity of Fe II emission with respect to other lines has not 
been well studied, especially in systematically for a large 
sample of objects (only a few m easurements of individua l 
sources have been published; e.g., IVeron-Cettv et al.ll2004l) . 
In fact, almost all broad quasar emission lines show blue- 
ward velocity shifts (Gaskell 1982; Carswell et al. 1991), with 
the exception of Mg II (no shift; Junkk arinen 1989) a nd H/ 3 
(redward shift in some studies; e.g., |McIntosh et al. !l999l) . 
Most previous studies simply assume d that Fe II has no 
shift with respective to [O III] (e.g., Boroson & Green||1992|^ 
Marziani et al. 1996; McLure & Jarvis 2002; Dietrich et al. 
2003^ .Greene & HQ,2005b: ,Kim et al.,.2006; ,Woo et aL,2006) 
and that it ha s the same line width as the broad component 
of up (e.g., ' Netzer & TrakhtenbrotI [20071; ISalviander et"ai] 
l2007h . The goal of this study is to test this assumption. 

The present paper presents the first detailed investigation 
of the velocity shift and width of optical Fe II emission^ in a 
large sample of quasars selected from the Sloan Digital Sky 
Survey (SDSS; York et ani200 0'). Our primary motivation is 
to determine the physical location and origin of the Fe II- 
emitting region. We describe selection of the sample in ^and 
spectral analysis in ^ We test the reliability of the measure- 
ments and the errors using Monte Carlo simulations 03.3l l, 
and then check how significantly our method improves the 
spectral fit and how our method affects the measurements of 
other emission-line parameters 03.4l i. Section|4]discusses the 
results we obtained, including the distribution of Fe II emis- 
sion shifts and widths, correlations with other parameters, and 
also an analysis of the composite spectra. The implications of 
our results are discussed in ^ with conclusing remarks given 
ing6] 

Throughout this work, we adopt the following cosmological 
parameters: Ho = 70 km s"' Mpc"', fi,,, = 0.3, and JIa = 0.7 
(ISpergel et aT]|2007h . 

2. SAMPLE SELECTION 

Our sample is selected from the SDSS Fifth Data Re- 
lease (DR5^_Ailelman-McCarthv et al. 2007) quasar cata- 
log ([Schneider et al.i [20071) . We choose objects with red- 
shift z < 0.8 to ensure that the [O III] emission line lies 
within the SDSS spectral coverage. Since the SDSS quasar 
sample is flux-limited and selected by broad-band colors 
(iRichards et al. 2002a), care must be exercised in using it 
to study the quasar lumin osity function ([Vanden Berk et alJ 
[20051; [Richards et al.[ [2006.) or its cosmological evolution. 
However, this sample is adequate for the scientific goals of 
this work. 

We impose a series of selection criteria to ensure reliable 
measurements of Fe II emission. (1) We require a signal- 
to-noise ratio (S/N) >10 in the wavelength range 4430-5550 

^ Unless otherwise noted, the Fe II emission in this paper refers to the 
optical band. 



A, covering H/3, [O III], and the most prominent features of 
optical Fe II emission. (2) We remove sources that have re- 
duced > 4 in the continuum decomposition 03.1.3I ). These 
sources cannot be fitted well by the present continuum model 
(Eq. ([T]) in ^3.11 ). (3) We remove sources that have EWpe < 
25 A; this EW cut is determined by the simulations described 
in ^3.31 (4) We also remove sources with H/3bc FWHM er- 
rors >10% and [O III] A5007 peak velocity shift errors >100 
km s"' . In total, the final sample contains 4037 sources, which 
is roughly 30% of all quasars in DR5 with redshift z < 0.8. 

3. SPECTRAL FITTING 

The spectra of quasars from UV to optical wavelengths are 
remarkably similar to each other. T he quasar composite spec- 
trum (e.g., Vanden B erk et alJ[2001[) is characterized by a fea- 
tureless continuum and a plethora of broad and narrow emis- 
sion lines. In the luminosity range of interest to us here, very 
little, if any, starlight is observed, so we can ignore the host 
galaxy contribution to the spectrum. The validity of this sim- 
plification can be tested from our spectral fitting — nearly all 
of the sources in the sample can be well fitted without a host 
galaxy component. 

The procedure of our spectral fitting algorithm is as follows. 
We begin by deredshifting the spectrum after correcting for 
Galactic extinction. Then, the continuum is decomposed into 
three components: (1) a single power law, (2) Balmer contin- 
uum emission (supplemented with high-order Balmer emis- 
sion lines), and (3) a pseudo-continuum due to blended Fe 
emission. We subtract this continuum model to obtain a pure 
emission-line spectrum, which is then fitted to derive parame- 
ters for the emission lines. The measurement of Fe emission is 
strongly affected by the uncertainty in the determination of the 
continuum level. Since there are almost no "pure" continuum 
windows, a simultaneous fit should be performed to decouple 
the Fe emission from the featureless continuum, rather than 
fitting the two independently. [Tsuzuki et al. (2006) adopt an 
alternative approach in which they use a theoretical model to 
estimate the flux fraction of the emission lines in the contin- 
uum windows.] The foUowing subsections will describe each 
step in detail. 

3.1. Continuum Decomposition 

We use the TJy-dependent Galactic extinction law given by 
'Cardefli et al. ( 1989) and assume Rv = 3.l. Eqs. (2a) and (2b) 
in .Cardelli et al.. (.1989.) are used for the infrared band, bu t 
Eqs. (3a) and (3b) are replaced by those in [0'Donnel ([1994 
for the optical band. We adopt the Galactic extinction in the u 
band listed in the SDSS quasar catalog ( Schneider et al. 2007^ 
and change it to the V band using the relation Ay = Auj 1 .579 
(ISchlegel etalJ[T998h . Then, we deredshift the extinction- 
corrected spectrum using the redshift provided by the SDSS 
pipeline. This value of redshift is later refined using the line 
centroid of [O III] A5007 measured after continuum decom- 
position (see 33.5l l. 

We model the continuum of the spectrum after Galactic 
extinction and redshift correction (hereinafter the corrected 
spectrum) using three components: 

+FP=(FFe,FWHMFe,yFe). (1) 

In total there are seven free parameters. The first term is the 
featureless power law 

^'^^ = ^^"'"(51^) ' 
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where F5100 is the flux density at 5100 A and a is the spectral 
index. The second and third terms denote the Balmer contin- 
uum and Fe emission, respectively, which are described be- 
low. For the Fe emission term, Fpe is the flux and Vpe is the 
shift velocity of Fe II. 

3.1.1. Balmer Continuum and High-order Balmer Lines 

Following iGrandil (Il982l) and iDietrich et all (l2002l) . flie 
Balmer continuum produced by a partially optically thick 
cloud with a uniform temperature can be expressed by 

Ff^ = FBEBA(7;)(l-e""^) (3) 

for wavelengths shortward of the Balmer edge (Abe = 3646 
A). B\{Te) is the Planck function at an electron temperature 
Te, FsE is a normalization coefficient for the flux at Abe, and 
Tx is the optical depth at A expressed by 

^^ = ^'^^(a^)' 

where tbe is the optical depth at the Balmer e dge. There are 
two free parameters, Fbe and tbe- Following IDietrich et alj 
( I2OO2.) , we assume T, = 15,000 K. 

At wavelengths A > Abe, blended higher-order Balmer lines 
give a smoot h rise i n the spectrum from ~4000 A to the 
Balmer edge dWills et al. 1985). We treat th e higher-order 
Balmer lines in a manner similar to that in IDietrich et al.l 
( |2003) . with some modifications. To determine the relative 
strengths of the transi tions with 1 <n < 50, we use the line 
emissivities given by IStorev & Hummer! (1 19951) for Case B, 
Te = 15,000 K, and = 10^ cm'l We normahze the flux of 
the higher-order Balmer lines to the flux of the Balmer contin- 
uum at the edge using the results in lWills eTal, (1985) . (This 
implies that the higher-order Balmer lines also depend on Fbe 
and Tbe ) In order to smooth the rise to the Balmer edge, we 
assume that each line has a Gaussian profile with FWHM = 
8000 km s~' . In practice, none of the assumptions concern- 
ing the higher-order Balmer lines actually impact our results 
because our fitting windows do not include this region (see 
EES- 

3.1.2. Fe Emission 

iPhillipsI (Il977h first introduced the template-fitting method 
to treat Fe emission in AGNs, using the Fe spectrum of the 
narrow-line Seyfert 1 (NLSl) galaxy I Zw 1 (z = 0.061) to 
construct an Fe II template. In most applications of this 
method, the amount of velocity broadening applied to the tem- 
plate during spectral fitting is either solved as a free parameter 
or is fixed (usually to the FWHM of broad H/3). But the Fe 
template itself is not allowed to shift in velocity. We follow 
essentially the same template-fitting, but we explicitly allow 
the width and shift of Fe II to be free parameters. We now 
describe the details of our algorithm. 

Considering the redshift range of our sample, we need both 
the UV and optical Fe template for I Zw 1 . In the UV band, we 
adopt the Fe template of Vestergaard & Wilkes (2001). Note 
that their template is set to zero around the Mg II line, which 
is unphysical. However, the Mg II line is not the main fo- 
cus of the present work, and for the purposes of this work 
we do not concern ourselves with this complication. In the 
optical, apart from the w idely used template constructed by 
iBoroson & GreenI (Il992h . some others are also available. For 
example, a template of any width c an be constructed f r om th e 
list of Fe Unes for I Zw 1 given bv lVeron-Cettv et all (|2004 . 



We compared the se two different template s and in the end 
chose the one from lBoroson & GreenI (1 19921) (kindly provided 
by T. A. Boroson) because it gives smaller reduced x^- How- 
ever, we have verified that the shift of Fe emission, one of the 
main goals of this work, is actually not sensitive to the choice 
of template. The velocity shifts measured using either of the 
two template are consistent with each other 

We combine the UV and optical templates to form a single 
template (with a gap from 3100 A to 3700 A, which has no 
data) and convolve it with a Gaussian function: 

Fl" = * G(Feonv,FWHMeo„v,Konv), (5) 

where /rizwi jg jj^g j 2w 1 Fe template, G is a Gaussian func- 
tion with flux Fconv, width FWHMconv, and peak velocity shift 
Vconv The convolution is done in logarithmic wavelength 
space because li(lnA) = dX/X = dv/c. The parameters in Eq. 
([T]) can be calculated as follows. The flux of the Fe emission, 
Ffb, is equal to F^om multiplied by the flux of the template. 
The shift of the Fe spectrum, Vfe, is simply Vconv Finally, the 
FWHM of the Fe lines can be expressed as 

FWHMpe = ^FWHM^^^i +FWHM2„„, . (6) 

In the above algorithm, we assume that the Fe emission in 
the UV and optical have the same width and velocity shift, 
and that the ratio of UV Fe flux to optical Fe flux is fixed to 
that of I Zw 1 . These assumptions help to reduce the number 
of free parameters, and seem appropriate given the S/N of the 
present sample. If w e split the UV/optical template following 
IVerner et al.l (|2004|) to three major wavelength bands — UV 
(2000-3000 A), small blue bump (3000-3500 A), and optical 
(4000-6000 A) — the final results will be determined mainly 
by the optical Fe emission. The reasons are as follows. First, 
the small b lue bum p Fe emission is outside of our fitting win- 
dows (see ^3. 1.3b . Second, the spectra of most of the quasars 
in the sample either do not cover or only cover a very narrow 
segment of the UV Fe emission wavelength range. Third, the 
S/N of the UV band is lower than the optical band. Thus, our 
measurements mainly trace the optical Fe emission and are 
not very sensitive to UV Fe emission. We consider each of 
the three assumptions in turn. 

The ratio of UV to optical Fe flux has been investigated 
by many authors (e .g., ISigut et al.ll2004t IVerneret all 120041: 
iBaldwin et al.l 12004 and references therein). This ratio de- 
pends on the physical parameters of the clouds, such as the 
hydrogen density nu, the hydrogen-ionizing flux <I>h, the ve- 
locity of turbulence, and so forth. However, for physical 
conditions ty pical of quasars (wh « 10" cm"^, ~ 3 x 
10^° cm~2 s~': lFerland et al.iri992l) . models of Fe emission in- 
dicate that there are large regions of parameter space where 
the ra tio is roughly constant (Figs. 3 and 4. in Verner et aT| 
l2004l) . This suggests that adopting a single ratio of UV to op- 
tical Fe flux (fixed to that of I Zw 1) should be a reasonably 
good approximation. 

The width of the UV Fe lines is also not necessarily equal 
to that of the optical Fe lines. Many authors fix the width of 
the optical Fe lines to the width of broad H/ 3bc, while in the 
UV th e width is fixed to that of Mg II (e.g.. ISalviander et al.l 
l2007h . This procedure implicitly assumes that the optical Fe 
lines and H/3bc originate from the same region, and simi- 
larly that the UV Fe lines follow Mg II. Empirically, how- 
ever, the width of H/3r c is consistent with that of Mg II 
dMcLure & Jarvisll2002l): the d ifference is only 0.05 dex on 
average dSalviander et al.l2007l) . Our results in ^4.2l also show 
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this. Thus, for the present purposes, it is safe to assume one 
single value for the width of the optical and UV Fe lines. We 
leave the width as a fr ee parameter, as has been d one in many 
previous studies (e.g., iMcLure & J arvis 2002; Di etrich et alj 
l20QllKimet"an '2006Y 

The velocity shift of the UV Fe lines^ could be different 
from that of the optical lines if they arise from different re- 
gions. We use only the optical Fe template and fit the con- 
tinuum only in the optical band for testing. The resulting Fe 
shifts change little, demonstrating that our measurements are 
mainly determined by the optical Fe lines. Thus, for our goal 
of studying the optical Fe lines, assuming one single shift for 
the UV and optical Fe lines is an adequate approximation. 

It should be noted here that because of the two reasons men- 
tioned above (low S/N around UV Fe lines and incomplete 
wavelength near), the present paper cannot conclude whether 
the UV and optical Fe II have a common origin. 

3.1.3. Multicomponent Fit 

The continuum model described by Eq. ([U is fitted by min- 
imizing the quantity 

\ 2 



where cr, is the error of the d ata set (;c,,y,). W e adopt the 
Levenberg-Marquardt method jPress et al.ll 19921 chap. 15.5) 
to solve Eq. ([T]i, which is nonlinear. We also use it in fitting 
the emission lines ( ^3.2l l. The fitting is performed in the fol- 
lowing windows: 2470-2625, 2675-2755, 2855-3010, 3625- 
3645, 4170^260, 4430-4770, 5080-5550, 6050-6200, and 
6890 -7010 A. These windows are devoid of strong emission 
fines dVestergaard & Wilkes 2001; Kim et al.l 2006). The win- 
dow 3625-3645 A is used to constrain the Balmer continuum 
emission, because in this region there is no strong Fe emission 
jWills et al.lll985h . The reduced distribution has a median 
value of = 1.365. 

Two examples of continuum decomposition are shown in 
Figures[T]and|2] The top panel shows the corrected spectrum. 
The spectrum in the fitting window is plotted in green. Each 
component is plotted in blue and the summed continuum in 
red. The middle panel shows the residual spectrum, which 
is the pure emission-line spectrum for the next step 03.2l i. 
The model fits the spectrum very well except in the region ~ 
3100-3700 A, where no Fe template is available. In fact, the 
total flux of the residual spectrum in this region strongly cor- 
relates with the Fe flux fpe. This is consistent with the "small 
blue bump" being produced by Fe lines and the Balmer con- 
tinuum (IWills et al.lll985l) . We subtract the power law and the 
Balmer continuum from the corrected spectrum and show the 
enlarged resultant spectrum (Fe-only spectrum) in the wave- 
length range 4100-5600 A in the bottom panel. The Fe model 
is plotted in red. 

From the first example on SDSS Jl 15507.61H-520129.6, 
which has narrow Fe lines, the Fe model not only agrees with 
the Fe-only spectrum in our fitting windows (in green), but it 
also fits the two strong Fe II lines at 4924 and 5018 A very 
well, even though these two lines are not in the fitting win- 
dows. Two blue dashed lines in the middle panel mark the 

^ In fact, the velocity shift of each individual UV Fe Une may be different. 
IVestergaard & Wilkes! j2001l) measured the velocity of each UV Fe line and 
showed that their shifts can vary by as much as ~ 100 km s"' . But consider- 
ing the resolution of our spectra, assuming a single velocity shift is safe, both 
for the UV and optical. 



positions of H/3 A4861 and [O III] A5007. The blue dashed 
line in the bottom panel marks the position of Fe II A4924 
line with zero velocity shift. The shift of Fe emission in this 
source can be seen clearly from the position of the Fe II A4924 
line; the velocity shift, as measured from the model and with 
respect to [O III], is 459±16 km s~'. For the second source, 
SDSS J112611.63H-425246.4, Fe II A4924 cannot be distin- 
guished from [O III] A4959 but the shift can be seen from the 
Fe II model. In this case, it is 1691±119km s"'. 

3.2. Emission-line Fitting 

After subtracting the continuum, we measure the H/3 and 
[O III] emission lines from the pure emission-line spectrum. 
We use multiple components to fit the emission lines over the 
wavelength range 4770-5080 A. The narrow H/3 component, 
[O III] A4959, and [O III] A5007 are modeled using three 
Gaussian. The [O III] AA4959, 5007 fines are forced to have 
the same FWHM and no relative wavelength shift, and their 
intensity ratio fixed to the theoretical value of 3.0. The nar- 
row H/3 component (H/3nc) is forced to have the same FWHM 
as [O III] A5007, a shift of up to 600 km s"' relative to [O 
III] A5007, and an intensity constrain ed to lie between 1/20 
and 1/3 of that of [O III] A5007 (e.g. lVeilleux & Osterbrocfl 
[19871 iMcGnreFal] 120081) . If necessary, we add another two 
Gaussian components for [O III] AA4959, 5007 to match their 
wings, and a corresponding Gaussian is added to H/3nc to 
ensure that H/3nc a nd [ O III] have the same profile. Fol- 
lowing |siivian^retaL| (l2003) and lMcGill et all (2008), the 
broa d H/3 component is modeled us ing a Gauss-Hermite func- 
tion (Ivan der Marel & Franxl[T993h . whose best fit yields the 
FWHM, peak velocity shift, and the square root of the second 
moment (cth/j). As illustration. Figure [3] shows the emission- 
line fitting for the two sources in Figures [T] and |2] The top 
panel shows the pure emission-line spectrum. The multiple 
components are in blue and the sum of them is in red; the 
bottom panel is the residuals. 

We also measure the [O ll] and Mg II emission lines when 
they are available. Since the continuum around [O ll] is 
not well determined (we make two assumptions in the model 
for the higher-order Balmer lines; see the last paragraph in 
^3.1.1b . we simply fit [O ll] above a loc ally defined contin- 
uum with a single Gaussian (iGreene & H o 2005a). As many 
sources have weak [O ll] emission, we adopt as detection 
criterion that the line must have an amplitude larger than 3 
times the standard deviation of the local continuum. The Mg 
II AA2796, 2803 doublet is fitted using two Gauss-Hermite 
functions; they have the same paramete rs except that the in - 
tensity ratio between them is fixed to 2 ( Baldwin et al.l ll996'). 
The FWHMs and velocity shifts of Mg II are calculated from 
a single Gauss-Hermite function. 

3.3. Tests of the Continuum Decomposition 

The template-fitting method for measuring Fe emission and 
the Levenberg-Marquardt algorithm for nonlinear fitting are 
widely used as almost "standard" approaches. We add a new 
parameter Vpe and set FWHMpe free in our fitting. Since the 
S/N of the majority of the sources in our sample are low 
(about 50% have S/N < 15), it is necessary to test the reli- 
ability of the Fe emission measuremen t. We perfo rm a suite 
of simulations similar to those done in iGreene & Ho ( 2006). 
Using a Monte Carlo method to generate artificial spectra, 
we measure the Fe emission of these spectra using the same 
method as that used for the observed spectra. Differences be- 
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Fig. 1. — Example of continuum decomposition for SDSS J115507. 61+520129. 6, which has nan'ow Fe lines. Its vpe is 459±16 km s '. Note that our Fe 
model fits the two strong Fe II lines at 4924 and 5018 A very well, even though these two lines are not in the fitting window. The top panel gives the spectrum 
after Galactic extinction and redshift correction. The spectrum in the fitting window is plotted in green; each component is plotted in blue; and the summed 
continuum is plotted in red. The middle panel shows the residual, pure emission-line spectrum. Two blue dashed lines mark the positions of H/3 A4861 and [O 
III] A5007. The bottom panel shows the spectrum after subtracting the power law and the Balmer continuum in the wavelength range 4100-5600 A. The red 
spectrum is our Fe model. The blue dashed line in the bottom panel is the position of the peak of Fe II A4924 with zero velocity shift. 



tween input and output parameters can then be compared to 
evaluate potential systematic errors and biases. 

We build the simulated continuum spectrum as a linear 
combination of a single power law and Fe emission expressed 
by Eq. (|5]). We generate a realistic noise pattern for the spectra 
using a real error array taken from SDSS observations, scal- 
ing it by a multiplicative factor to match the desired S/N of the 
simulation. For each pixel of the simulated spectrum, a Gaus- 
sian random deviation is added. Bad pixels have large devi- 
ates statistically. We use a mask array from an actual FITS file 
to locate the masked pixels; some of the pixels that have large 
errors are masked by the SDSS pipeline, but not all. This pro- 



cedure ensures that the simulated spectra have a realistic noise 
level and noise pattern. 

There are three main factors that can affect the measure- 
ments of Vpe and FWHMpe: (1) the S/N of the spectrum; 
(2) the strength (EW) of the iron emission; and (3) the width 
(FWHM) of the iron lines. Our simulations demonstrate that 
the input value of Vpe in the simulated spectrum has a very 
minimal effect on the systematic bias of the measured output 
values of Vpe or FWHMpe, and so we neglect it from further 
consideration. We calculate EWpe = F(Fe II A4570)/f5ioo 
where, in the present paper, F(Fe II A4570) is the flux of the 
Fe II emission between 4434 and 4684 A. Most of the quasars 
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Fig. 2. — Same as Fig. Q] but for SDSS Jl 1261 1.63+425246.4, which has broad Fe hnes. Note the large departure of the modeled Fe II A4924 line from its 
rest-frame wavelength (blue dashed line in the bottom panel). I'pe is 1691±1 19 km s"' . 



in the sample have a EWpe between 15 and 75 A and FWHMpe 
between 1000 and 5000 km s~' . Accordingly, we test four val- 
ues of EWpe (15, 25, 50, 75 A) and vary FWHMpe from 1000 
to 5000 km s"', in steps of 500 km s'K We set S/N = 10, 
which is the lower limit of the S/N in the sample. Vpe is fixed 
at zero so that we can examine whether the measured shifts 
are real or spurious. 

For each pair of values for EWpe and FWHMpe, we gener- 
ate 100 spectra. We fit the continuum and then calculate the 
quantity 

""■"FWHMin' 

where Vi„ and Vout are input and output values of Vpe, respec- 
tively, and FWHMin is the input FWHMpe- Figure IIJ; shows 
S^^^ as a function of FWHMpe- The results are plotted using 



different colors and line styles for different values of EWpe- 
For each EWpe, there are three lines: the middle line repre- 
sents the mean S^^^, and the other two hnes above and below 
are one standard deviations (aj^j^) above and below the mean. 
From the diagram, we can see that the average value of (5^^ 
is always close to 0. This demonstrates that the measured 
value of Vpe shows no systematic redshift or blueshift with in- 
creasing FWHMpe and EWpe. Obviously, cr^j^ increases with 
decreasing EWpe; measuring Fe II emission is very uncertain 
when EWpe is small. Based on the results of these tests, we 
decided to exclude from the sample sources with EWpe < 25 
A. It should be noted that the S/N in this particular simulation 
is set at the lower limit of our sample, so the uncertainties 
shown in Figure |4] should be considered upper limits. Fig- 
ure I4J7 shows the difference between the error given by our 
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Rest Wovelength (A) Rest Wovelength (A) 

Fig. 3. — Fitting of the emission lines for the two quasars in Figs.[T]and[2] The top panel shows the pure emission-Hne spectrum. Each component is in blue, 
and the sum of them is in red. The bottom panel shows the residuals. 



code and that given by the simulation, a^^ is the average value 
of the errors given by the code in the 100 trials divided by 
FWHMi„. 

FiguresSJ: andHJ:/ show the systematic error in the measure- 
ment of FWHMpe- We define 

cFWHM _ FWHMput ^ 

^™ FWHMin ■ ^ ' 

The standard deviation of f^!^^'^ is denoted by cr^™. 



^FWHM 



and 

is the average value of the relative errors of FWHMpe 
(defined as the error of the FWHMpe divided by FWHMin) 
given by the code. 

From these simulations, we conclude: (1) the mean values 
°f '^sim '^sim™ cluster around and exhibit no trend as 
a function of FWHMpe or EWpe; (2) the error given by our 
code is consistent with that given by the simulations (except 
perhaps for the EWpe = 15 A bin). These results show that our 
measurements are reliable and robust. 

3.4. Comparison with Conventional Fe Template-fitting 
Methods 

There are two conventional methods for Fe fitting. 
Both have been widely used and effective. The first 
(model 1) assumes that the Fe lines have no velocity 
shift but that their width can be different with that of 
H/3bc (e.g., Boroson & Green 1992; Marziani et al. 1996; 
McLure & Jarvi s 2002; Dietr ich et all 120031: IGreene & Hoi 
2005bt iKim'eTal., .2006; Woo et al.1 l2006h . Alternatively 
(model 2), one assumes that the Fe Unes have no 
shift and th at they have the same width as H/3b c (e-g-, 
iNetzer & T rakhtenbrot 2007; Salviande Tetal] l2007t). Two 
obvious questions arise. Does our continuum model (Eq. 
([TJ) significantly improve the fit? And second, how does our 
continuum decomposition affect the emission-line measure- 
ments (e.g., H/3 and [O III]) and the physical parameters sub- 
sequently derived from them (e.g., central BH mass and the 
Eddington ratio)? 

We compare the results derived from our approach with the 
two standard methods described above. Following iLuptonI 



CHI Chapter 12.1), we use the F-test^ to calculate how sig- 
nificantly our model improves the fit for each source. Com- 
paring with model 1, 71% of the sources are better fit by our 
model at a significance of >95.45%, and 58% of the sources 
are better fit at a significance level of >99.73%. With respect 
to model 2, the corresponding improvement can be seen in 
89% of the sources at a significance of >95.45% and in 80% 
of the sources at a significance of >99.73%. On average, our 
model decreases the reduced by 0.039 and 0.094 compared 
with models 1 and 2, respectively. We conclude that our ap- 
proach of allowing Vpe and FWHMpe to be free parameters 
significantly improves the fit in most objects. 

Next, we evaluate the actual impact that the different meth- 
ods have on measured and derived physical quantities. We 
measure the H/3 bc an d [O III] emission lines, derive Mbh and 
Lboi/iEdd (see 34.41 ) for each model, and then calculate the 
relative differences of the parameters between our model and 
the two fiducial standard models. As summarized in Table[Tl 
the differences in line luminosities and line widths for H/3bc 
and [O III], Mbh, and Lboi/i'Edd are all less than 5%, while 
changes in velocity shifts are also no more than 50 km s"' . 
The only exception is for sources with very weak [O III] lines. 
In this regime, the measurement of [O III] can be strongly 
affected by Fe II A4924 and Fe II A5018, and the effect on 
L|-Q jjjj and FWHM|^q jjjj in Table[T]is large (a few tens of 
percent). This exercise demonstrates that, for most applica- 
tions, the choice of method for Fe template fitting is in prac- 
tice unimportant — unless the main scientific objective is to 
actually study the Fe II emission itself. 

3.5. Redshifts 

The narrow emission lines are commonly used to obtain the 
systemic redshift. The [O III] A5007 line is the strongest nar- 
row line for mo st quasars, so it is most often used. However, 
iBorosonI (l2005h showed that [O III] can be blueshifted with 

^ Strictly speaking, this test is valid only for models that use linear fit- 
ting. But it has be en empirically u sed for nonlinear models and seems to be 
effective (e.g., see lHao et al.l2003) . 
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FWHMpg (km s ^) FWHMp^ (km s '') 

Fig. 4. — (a) Input-output simulations of vpc measurements. We define 5^^^ = (Vin — Voiii)/FWHMi„. For eacli EWpc, the middle line represents the mean S^^^, 
and the other two lines above and below are one standard deviation above and below the mean. The S/N in the simulation is 10, which is the lower hmit of our 
sample, so the uncertainties shown here should be considered upper limits, (b) Test of the error bai's given by the code, is the average value of the errors given 
by the code in 100 trials, divided by FWHMin, and ct^^^ is the standard deviation of 5^^. (c) Input-output simulations of FWUMpe measurements. We define 
^FWHM ^ (FWHMou, -FWHMin)/FWHMi„. (d) Test of the eiTor bars by the code, o-™"^^ is the average value of the relative errors of FWHMpc given by the 
code in 100 trials, and o"™'^'^ is the standard deviation of fSjj™'^. See text for details of the simulation. 



TABLE 1 

Effect of Fe Template-fitting Method on Other Parameters 



Model 
(1) 


(2) 


FWHM(H/3bc) 
(3) 


V(H/3bc) 
(4) 


^[Olll] 
(5) 


FWHM[oi„] 
(6) 


^[0 III] 
(V) 


Mbh 
(8) 


^bol/^-Edd 

(9) 


fix Vfc, free FWHMpe 
fix both Vfc and FWHMpc 


-1.75%(2.99%) 
-1.46%(3.14%) 


-1.74%(3.54%) 
-0.39%(4.10%) 


-27.2(49.2) 
-23.2(47.9) 


-2.15%(5.01%) 
-0.48%(5.58%) 


-1.02%(9.51%) 
-0.31%(10.22%) 


2.76(35.6) 
2.95(32.0) 


-0.07%(0.49%) 
-0.01%(0.46%) 


2.23%(11.2%) 
0.63%(10.5%) 



Note. — Changes in emission-line parameters and derived physical parameters due to using different Fe template-fitting models. Col. (2): luminosity of H/3bc- 
Col. (3): FWHM of H/3bc- Col. (4): H/3bc velocity shift. Col. (5): luminosity of [O III]. Col. (6): FWHM of [O III]. Col. (7): [O III] velocity shift. Col. (8): mass 
of the central BH. Col. (9): Eddington ratio of the central BH (see j|4.4l for details of how to derive Mbh and ibol/^Edd)- Each column shows the relative changes in 
percentage, except for Col. (4) and (7), which show the difference in absolute velocity shift in km s"' . The number in parenthesis is the standard deviation. 
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[0 lll]-[0 II] Velocity Shift (km s~^) 

Fig. 5. — Distribution of [O III] velocity shifts with respect to [O II]. 
Positive velocity indicates a redshift. The dotted line marks the position of 
zero velocity shift, while the dashed line marks the median of the distribution. 



respect to the low-ionization forbidden lines, which provide 
a better rest-frame. In our sample, 2265 sources 50%) 
have detectable [O ll] A3727 emission. Figure |5] shows the 
distribution of velocity shifts between [O III] and [O ll]; pos- 
itive velocit y indic a tes a redward shift. Consistent with the 
results of Boroson (2005), the majority of the sources have 
blueshifted [O III]; the median blueshift is -47 km s"'. Be- 
cause the [O III] blueshifts are much smaller tha n the velocity 
shifts seen in Fe II and in other broad lines (see ^4.11) and only 
half of our sample is detected in [O ll], we still use [O III] to 
define the rest-frame. Thus, we define the Fe II velocity shift 
by VFe=VFe-V[Q jjj] , whcrc Vj-Q jjj] is the velocity shift of the 
core component of the [O III] A5007 line measured in ^3.21 

3.6. Errors 

We calculate the error of vpe from the fitting of the contin- 
uum and the fitting of the [O III] line. Figure |6]5( shows the 
distribution of the vpe errors. Most velocity shifts have an er- 
ror < 200 km s"', and the median of the distribution is 116 
km s~' . We also plot the relative error on FWHMpe in Figure 
lit'. The typical value is about 10% to 20%, and the median is 
12.0%. 

4. RESULTS 

The full catalog of the measurements used in the analysis 
below is available electronically. Table |2] describes the con- 
tents and the formats of each column. All the emission line 
velocity shifts are with respect to [O III] A5007 and a positive 
value indicates a redward shift. -9999 in the columns of Mg 
II (Col. 21—26) and [O ll] (Col. 35 — 40) measurements in- 
dicates the line lies out of the SDSS spectral coverage or too 
weak to be detected (see ^3.2l for details). The details of how 
to deriving the radio and X-ray properties (Col. 19 & 20), and 
the Mbh and Eddington ratio Lboi/i'Edd (Col. 42 & 43), are 
described in ^4.5l and ^4.4h 'espectivelv. 

4.1. Fe Emission Shifts 



From the distribution of vpe (Fig. [TJ?), we can clearly see 
that most quasars exhibit Fe II emission that is redshifted 
with respect to the systemic velocity of the narrow-line re- 
gion (defined by [O III]). Only 481 out of 4037 quasars have 
blueshifts. The median shift is vpe = +407 km s"' (the verti- 
cal dashed line in the figure), with a maximum value of vpe 
« 2000 km s"'. Considering that the typical error on vpe is 
only <200 km s"' 03.6I ). the vast majority of the values in 
the redshifted tail of the vpe distribution must be real. For 
comparison, we also calculate the distribution of Fe II veloc- 
ity shift with respect to [O ll] (Fig. IT!?). The results obtained 
by using [O ll] as reference instead of [O III] are very similar 
and lead to the same conclusion. 

To rule out the possibility that the excess Fe II redshifts 
arise from artifacts due to poor data quality or fitting errors, 
we examined a subset of data using the much stricter selection 
criteria that S/N > 15 and 1-2. These 309 quasars have 
the best data quality and the most reliable continuum fitting. 
The shape of the distribution of vpe this subsample (Fig. |7}:) is 
quite similar to that of our whole sample. The median velocity 
shift is 414 km s"'. We conclude on the basis of this test, as 
well as the Monte Carlo simulations described in ^3.31 that 
the Fe II redshifts are robust and reliable. 

Finally, we illustrate that the velocity shifts found for Fe II 
really do imply radial motions of the Fe ll-emitting region, 
and not others. Figures|7j/ and|7^ show the velocity shifts for 
H/3bc [v'(H/3bc)] and Mg II. The distribution of v(H/3bc) is 
almost symmetrical around 0, with a median value of only 
65 km s"' and a maximum value of about ±1000 km s"'. 
This result is consistent with previous studies (e.g.. Fig. 3 of 
Sulentic et al. 2000a; Fig. 2 of Baskin & Laor 2005; Fig. 6 of 
Isiiang et al. 2007). Most recently. iBonning et all (1200^ also 
studied the velocity shift of H/9bc using SDSS quasars, and 
our result is very similar to theirs; they fit their distribution of 
velocity shifts (their Fig. 1) using a Gaussian profile with a 
peak velocity of 100 km s"' . Our distribution of Mg II shifts is 
indistinguishable from that of H/9bc- It is symmetrical around 
0, has a median value of 1 13 km s"' and a maximum value of 
about ±100 km s~\ and it is consi sten t with those given i n, 
for example. [Richards et alJ (l2002bh and lShang et al.r(l2007h . 

In order to compare Fe II directly with H/9bc, we also plot 
the distribution of the Fe II velocity shift with respect to H/3bc, 
as shown in Figure |7f. Again, most objects have redward 
shifts in Fe II. 

4.2. Fe Emission Widths 

Figures [H]^ and [8t> show the distribution of FWHMpe and 
FWHM(H/?bc), respectively. FWHMpe has an artifical lower 
limit of 900 km s"' , which is bounded by the I Zw 1 Fe II 
template. The median value of FWHMpe, 2533 km s"', is 
- 0.74 of that of FWHM(H/3bc) (3445 km s'^). Almost all 
the sources have Fe II lines narrower than H/3bc, and the ma- 
jority have FWHMpe w 3/4 FWHM(H/3bc) (Fig. [Hh). This 
result is contrary to the prevailing notion that Fe II and H/3bc 
have similar pro files an d are emitted from the same region 
(e.g.. Boroson & Greenl lT992). but is cons istent with some 
more recent stud ies of a few objects (e.g.. iPopovic et all2007l ; 
iMatsuoka et"ani200&) . 

We find that FWHM(H/3) and FWHM(Mg ll) are 
well correlated and roughly equ al (Fig. [8]f), consis - 
tent wit h McLure & Jarvis (2002); ISalviander et all (l2007h . 
ISalviander e t al. (2007) find that FWHM(H/?) tends to be 
larger than FWHM(Mg ll) for FWHM(H^) > 4000 km s"', 
an effect they attribute to an extensive red wing on H/3. This 
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Error of the Fe Velocity Shift (km s ) Relative Error of the Fe FWHM 

Distribution of (a) error in tlie Fe II velocity siiift and (b) relative error in Fe II FWHM. The dashed line marks the median of the distribution. 



tendency can also be seen in our plot, and our composite spec- 
tra (Fig. [T3T l do show that H/3 tends to have a red asymmetry 
when broad. 

4.3. Correlations with Other Emission-line Parameters 



iBoroson & GreenI (|1992|) used principal components anal- 
ysis to study the correlations among various observed prop- 
erties of nearby quasars, and found that most of the variance 
in the optical spectra of quasars is connected with the inverse 
coiTelation between Fe II and [O III] strength. The soft X- 
ray photon index (Fsoft), the ratio of Fe II to H/3bc (^Fe), and 
the h ne width of H/3bc fFWHM(H^Bc)1 correlate with each 



other jBoller et al.ll996l;jWang et al.ll996HLaor et al.lll997al) . 
iSulentic et all (l2000albL |2007l) identified that FWHM(H/3bc), 
/?Fe, Fsoft, and the velocity shift at half maximum of the broad 
C IV line profile [c(i)] provide discrimination between dif- 
ferent AGN types. So, for simplicity, instead of using a for- 
mal principal components analysis to study the correlations 
between Fe II and other parameters, we just study the coiTela- 
tions with FWHM(H/3bc) and /Jpe- 

The left panel of Figure |9] shows the dependence of vpe 
on FWHM(H/3bc) for our sample. The crosses in the figure 
present the median values of the errors in both coordinates. 
Dividing the sample by a vertical line of FWHM(H/3bc) = 
3000 km s"' and a horizontal line of vpe = 800 km s"' , we find 
that almost all sources with FWHM(H/3bc) < 3000 km s"' 
have vpe < 800 km s"', and most sources with vpe > 800 
km s"' have FWHM(H/3bc) > 3000 km s"'. Sources with 
large vpe and low FWHM(H/3bc) are rare. There is also a 
large spread in vpe for sources with FWHM(H/?bc) > 3000 
km s"'. Note that naiTow-line Seyfert 1 galaxies (NLSls), de- 
fined by FWHM(H/3bc)< 2000 km s"' (l6sterbrock'& Pogg^ 
Il985h . almo st all have low vp^ in our sampl e. 

Following iNetzer & Trakhtenbroj (l2007l) . we define 7?Pe as 
the line luminosity ratio of Fe II and H/3bc, 



tween A4434 and A4684. The right panel of Figure |9] shows 
the correlation diagram for vpe versus R^^. The vertical line is 
7?pe = 1. Most sources with /?pe > 1 have small vpe- Sources 
with large R^^ and large vpe are rare. This correlation agrees 
with the correlation between /?pe and Lboi/^Edd shown in Fig- 
ure 5 of Netzer & Trakhtenbrot (2007), considering the strong 
inverse correlation between vpe and Lboi/i'Edd as we find in 
g4l below. 

The statistical connection between vpe and FWHM(H/3bc) 
or j^pe is s imilar to those shown in Figures 2 and 9 of 
IBoroson &~ Green ( 1992). This suggests that vpe can also pro- 
vide useful empirical discrimination between different types 
of AGNs from optical spectra. 

4.4. The Physical Driver ofv^e 

The Eddington ratio Lboi/^Edd is often suggested to be 
the main physical dri ve r of the spectral diversity in 



(ISulentic et al. 



iMarziani et al. 



ical dri ve r or the spectral diversity m quasars 
I2000allbl: iMarziani et al.ll200H lBorosonil2002l: 



12003 al) . and Mbh is argued to be an impor- 



tant determinant of radio-loudness (Laor' 20()0tlBorosonl2002l; 
McLure & Jarvis 2004, but see Ho 2002). This section inves- 
tigates attempts to determine which physical variable is the 
main driver of variations in vpe. 

Central black hole masses can be estimated from empiri- 
cal relations derived from reverberation mapping. We derived 
the BH mass and L boi/Lpdd using the relation calibrated by 
iMcGill etalJ (120081) . 



log 



M, 



BH 



Mr. 



= 7.383 + 2 log 
+0.69 log 



1000 km S-' 
Aisioo 



1044 



ergs s 



(11) 



Rfp = 



L(Fe II A4570) 



(10) 



where L(Fe II A4570) is the luminosity of Fe II emission be- 



(their Table 3; we use the factors for L5100./ and anp)- 
We est imate the bolometric luminosity using Lboi = 9AL5100 
dKaspi et al. 2000). 

FigurefTOb shows vpe as a function of Lboi/i'Edd for our sam- 
ple. The vertical dot-dashed line is log(Lboi/iEdd) = -0.8, 
and the horizontal dot-dashed line is vpe = 800 km s"'. We 
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Fig. 7. — (a) Distribution of Fe II velocity shifts witli respect to [O III], (h) Fe II velocity shifts with respect to [O II]. Panel (c) shows the distributions of 
vpe in the subsample that has stricter criteria S/N > 15 and < 1-2. (d) H/3bc velocity shifts with respect to [O III], (e) Mg II velocity shifts with respect to 
[O III] . (f) Distribution of Fe II velocity shifts with respect to H/3bc- Positive velocity indicates a redshift. The number in each panel is the median value of the 
respective distribution, which is also marked by the dashed line. The dotted line marks the position of zero velocity shift. 
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TABLE 2 

Spectrophotometric Measurements Table Format 



Column 


Format 


Units 


Label 


Description 


1 


al8 




SDSS Name 


SDSS DR5 Object Designation hhmmss.ss+ddmmss.s (J2000.0) 


2 


f8.4 




z 


Redshift determined using the peak of the [0 III] A5007 


3 


el 2.4 


g|*gg g 1 




T iiminr»cit\? r\f thf T?** IT f»inioc;ir\n V\f*hj/Af»n AA'2A ctnfi A.f\^f\ A 
LjLllllillUftil y yJi LllC 1 C li CIIUoolUlI UC-IWCC-II ' l ' l J'-t ttllU i\. 


4 


el 2.4 






Error in L(Fe II A4570) 


5 


fg 1 


km s"' 


A YT J. J-lVAJng 


Fe IT FWHM 

1. C 11 1. V¥ lJ.iVl 


6 


f8.1 






Error in FWHMpe 


7 
8 


fg 1 
f8.1 


km s~' 


Vpe 


Pf* TT vflnr'itv *;niTf' nil thf* fmiQQinn linf* vflnpitv Qnift*; nff* 

1. W 11 Vt^lU^lLy jlllll., U.11 Lllw WlllljjlVJll llllw V^lU^lLy JllllLS til t 

with respect to [0 III] A5007, and a positive value indicates a redward shift. 
Error in vpe 


9 


el2.4 


ergs s~' 


i(H/3Bc) 


H/3bc luminosity 


10 


el2.4 






Error in L(H/3bc) 


11 


f8.1 


km s~' 


FWHM(H/3bc) 


H/3bc FWHM 


12 


f8.1 






Error in FWHM(H/3bc) 


13 


f8.1 


km s~' 


v(H/3bc) 


H/3bc velocity shift 


14 


f8.1 






Error in v(H/3bc) 


ID 


eiz.4 


ergs s~' 


T 

i'5100 


Specific continuum luminosity at 5 100 A 


ID 


eiz.M- 






Error in /^sioo 


1 7 

1 / 


f7 'X 

L / .J 




cx 


Power law spectral index of the continuum 


18 


t7.3 






Error in cx 


19 


a2el0.3 


ergs s^^ Hz^^ 




Specific luminosity at 6 cm derived from the Peak flux density measured in FIRST, 
assuming a radio spectral index = -0.5 ; the upper limits are derived from 
the FIRST flux limits; -1.000 indicates not in FIRST survey area. 


20 


ell.3 


ergs s~' Hz~' 


i^keV 
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Fig. 8. — Distribution of (a) Fe 11 FWHM and (h) H/3bc FWHM. Tlie number in the panel is the median FWHM, which is also marked by the dashed line. 
FWHMpc has a lower limit of 900 km s"' (dotted line); this is artificial because the I Zw 1 Fe template we used has a width of 900 km s"' . (f ) Correlation 
between FWHMpe and FWHM(H/3bc)- The dashed diagonal line shows FWHMpc = FWHM(H/3bc), and the solid line is FWHMpe = 3/4 FWHM(H/3bc)- (d) 
Correlation between FWHMpe and FWHM(Mg II). The dashed diagonal denotes FWHM(H/3bc) = FWHM(Mg II). 

an v'Fe of about 150 km s~ ', th e scatter to the fitted line in- 
creases. As mentioned in 93.31 the input value of Vpe affects 
the measurements little; this means that the errors of vpe will 
not decrease with vpe- Thus, for sources with small vpe, the 
fractional error on vpe will be large. This causes the scat- 
ter described before. The squares show the mean values of 
log(Lboi/i'Edd) in bins of Avpe = 200 km s"'; the error bars in 
panel (a) show the standard deviations. The dashed line shows 
the fit to the binned values of log(Lboi/i'Edd): 



find that almost all sources with log(Lboi/i'Edd) > -0.8 have 
vpe < 800 km s"', and most sources with vpe > 800 km s"' 
have log(Lboi/i-Edd) < -0.8. There are very few sources with 
large vpe and lai-ge Lboi/ipdd- The eiTor in Lboi/^^Edd is roughly 
0.3 dex, which is dominated by the systematica l error in es- 
timati ng the BH mass using empirical relations (IMcGill et al.l 
l2008h . 

The logarithmic form of this diagram is plotted in Figure 
\Wb . excluding the 481 sources with negative vpe. We find a 
strong inverse correlation between vpe and Lboi/^Edd- The fit 
taken into account the uncertainties in both quantities yields 
the solid line: 

log ype = (1 .00 ± 0.05) - (1 .83 ± 0.05) log(Lboi/iEdd). (12) 

Pearson's correlation coefficient rp is -0.53, and the proba- 
bility P of a chance correlation < 1 x 10"^. Note that below 



log ype = (0.60 ±1.40)- (2.22 ± 0.99) log(Lboi/iEdd) ,(13) 

In this case, rp = -0.98 and f < 1 x 10"''. The above 
analysis indicates that vpe depends strongly on Lboi/i^Edd; 
vpe oc (iboi/i'Edd)^, with 7 « -2: the larger the Eddington 
ratio, the lower the velocity shift vpe. Plotting vpe as an func- 
tion of Mbh and Lboi (not shown) reveals that neither of these 




log(^boiAEdd) i°g(^boiAEdd) 

Fig. 10. — Correlation of Fe 11 velocity shift vs. Eddington ratio, with the Fe II velocity shift plotted (a) linearly and (b) logarithmically. See text explanation 
of the dot-dashed lines. The solid lines is the fit to data in Eq. )12t . The green squares represent the mean log(Lbol/iEdd) in bins of Avpe = 200 km s"' , and the 
error bars represent the standard deviations. The green dashed fine shows the fit to the binned data. 



two variables is as important as Lboi/ifidd- The Eddington ra- 
tio is the main physical driver for vpe- This result provides a 
strong constraint on theoretical models of the Fe II emission 
region. 

4.5. Correlations with Radio and X-ray Properties 

In an effort to understand the physical origin of the Fe II 
velocity shift, we examine whether vpe cor relates with radio 
and X-ray emission. Richards et all (l2002bl) conducted a sim- 
ilar investigation in their analysis of velocity shifts for the C 
IV line. Figure [TTI (upper panel) plots the fraction of radio- 
loud quasars in bins of different vpe. We define the radio- 



loudness parameter as = log(L6cm/iB), where Lecm and Lb 
are the obs erved luminosities at 6 cm and 4400 A. We use 
the FIRST (^Becker et al.l| I995h peak flux densities at 20 cm 
in Table 2 of Schneider et alj ( l2007h to calculate the radio 
luminosity, assuming a radio spectral index = -0.5. The 
optical luminosities are calculated from the power-law con- 
tinuum we fitted. We classify the sources with R> I as radio- 
loud. There are a total of 165 radio-loud quasars out of 3750 
sources in the present sample within the FIRST survey area. 
Note that the percentage of radio-loud quasars in our sam- 
ple is 165/3750 = 4.4%, slightly lower than the 6.7% found 
bv lMcLure~& JarvisI ( |2004 . The reason is probably that our 
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sample is biased toward quasars with high Lboi/iEdd (by our 
Fe II EW cut) , since decr eases with increasing Lboi/i'Edd 
(lHoll2002t IGreene et alj 120061) . The rightmost bin have the 
largest radio-loud fraction, in which there are 7 radio-loud 
quasars out of 85, a fraction ~2 times higher than average. 
The cumulative Poisson probability (P) of getting 7 or more 
objects out of 85 is 8.54 x 10"^ (Gehrels 1986) when the aver- 
age is 3.74 (165/3750 x 85), no more than a 3(t significance 
(1 -f < 99.73%). Thus there is no clear trend of radio-loud 
fraction with vpe- 

Next, we evaluate the fraction of X-ray-detected quasars 
as a function of vpe (middle pan el of Fig. fTTT l. We use 
the data from S chneider et aTl (2007), who provide the X-ray 
full-band count ra te from the ROSAT All-Sky S urvey Bright 
(IVoges etal.ll 19991) and Faint (IVoges et al.ll2000l) sources cat- 
alogs. No obvious trend is apparent. Finally, we test for pos- 
sible dependence of vpe on the optical-to-X-ray spectral index 
aox = -0.3838 log(L25{)o/^2 kev). where L2500 is the specific 
luminosity at 2500 A calculated using the power-law contin- 
uum we measured and L2 kev is the specific luminos ity at 2 
keV d erived from the ROSAT count rate using P IMMS (iMukail 
[T993h assuming a power-law model with photon index of 2. 
We see no obvious trend between aox and vpe either (bottom 
panel of Fig. fTTT l. 

4.6. Composite Spectra 

In order to get a visual impression of the correlations be- 
tween vpe and other emission-line properties, we create a set 
of five composite spectra by combining quasars in bins of 
different vpe. We divide our sources into five subsamples, 
covering the following velocity ranges: -250 to 250 km s"' 
(A, 1350 objects), 250 to 750 km s"' (B, 1362 objects), 750 
to 1250 km s-i (C, 590 objects), 1250 to 1750 km s"' (D, 
332 objects), and 1750 to 2250 km s"' (E, 180 objects). The 
composite spectra are generated following the procedure of 
IVanden Berk et a L (2001). The spectra of quasars in each sub- 
sample are deredshifted using the redshifts determined from 
[O III] and then normalized to unity average flux density over 



the rest wavelength interval 5090-5110 A. We generate the 
composites using the geometric mean, which is appropriate 
for quasars with power-law spectra because the geometric 
mean will result in a pow er law with the mean spectral index 
dVanden Be rk et al.ll2001h . 

Figure [12] shows the composite spectra of the five subsam- 
ples described above. We plot the spectra in different colors 
and shifted them slightly vertically for clarity. The spectra are 
arranged so that, from top to bottom, the velocity shifts of the 
Fe II emission increase. We find that FWHM(H/3bc) increases 
while the Fe II flux decreases, consistent with the correla- 
tions seen in 94.31 Inspecting the continuum slope, composite 
A shows a redder continuum than the rest (this can be seen 
most easily when comparing composites A with B). This sug- 
gests that sources with low vpe, which from our analysis are 
often accompanied by high Lboi/i^Edd and strong Fe II emis- 
sion, tend to have redder UV-optical continua. This pat tern is 
reminiscent of that seen bv lConstantin & Shieldsl ( l2003h . who 
found that the UV-optical continuum of NLSls, which have 
high Lboi/i'Edd and strong Fe II, is redder than that of regu- 
lar AGNs. This is unexpected from standard accretion disk 
models, since a higher Lboi/ipdd ge nerally produces a h otter 
disk and thus a bluer continuum (e.g.. lHubeny et alJ2000L and 
references therein). 

Note that it is difficult to see the shifts of the Fe II emission 
clearly in the composite spectra. For example, the peaks be- 
tween ^5150 and 5250 A appear to remain unshifted. The 
reason is as follows. On the one hand, since both the shifts 
and widths of the Fe II emission have a large scatter in each 
subsample, stacking the spectra effectively smooths the Fe II 
emission. On the other hand , there are many narrow e mission 
lines (e.g., line system N3 in lVeron-Cetty et ani2004l) that are 
weak in a single spectrum and have nearly no shift. Stacking 
the spectra enhances these narrow features. 

Another interesting phenomenon is that the wings of the H/3 
profile become progressively more (red) asymmetric when 
vpe increases (Fig. [T3T i. The inserted plot shows the differ- 
ence spectra between composites B through E, using com- 
posite A as reference, to emphasize the profile changes. 
Note the systematic migration of the red excess as vpe in- 
creases. We do find so me individual sources, similar to 
OQ 208 (iMarziani et all i 1993). whose redshifted H/3 com- 
ponent seems to be associated with Fe 11 emission ( 95.1b . 
The asymmetry in H/3 h as been studied by many authors. 
iBoroson & GreenI ([1992') found that there are H/3 red asym- 
metries at small iJpe. Marziani et al. (1996) found that radio- 
loud AGNs show predominantly redshifted and red asymmet- 
ric H/3 profiles: the larger the shift, th e broader the line. Re- 
cently, iNetzer & Trakhtenbroj (l2007h investigated the frac- 
tional luminosity of the red part of the H/3 line and found that 
it has a tendency to increase with decreasing Fe ll/H/3 (see 
their Table 3). Considering that vpe inversely correlates with 
7?pe, as shown in the right panel of Figure|9] the finding here is 
consistent with theirs. One possible interpretation, discussed 
in 95.21 is that the H/3 excess emerges from the same region 
that produces Fe II. A systematic study of this issue will be 
carried out in a future paper. 

5. DISCUSSION 

5.1. Redshifted H[3 Component: OQ 208-like Quasars 

In our sample, we find a class of sources whose H/3 pro- 
file can be fitted well including an additional, substantially 
redshifted Gaussian component. Usually this additional com- 
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Fig. 12. — Composite spectra of quasars in bins of different vpe. See text for details on generating the composite spectra. We shifted the spectra vertically for 
claiity. Major emission Hnes are indicated; see lVanden Berk et alj 4200TI) for a complete identification of the emission lines. 



ponent has a velocity width intermediate between that of the 
broad and narrow components. We call this the intermedi- 
ate component. An interesting and possibly highly significant 
fact is that the width and velocity shift of this additional H/3 
component is consistent with those of the Fe II emission. A 
prototype of this kind of sources is OQ 208 (M rk 668 ), which 
was first studied in detail by Marziani et jil] ( Il993h . They 
pointed out that in OQ 208 the Fe II hnes at 4924 and 5018 A 
have the same peak displacement as the red peak of H(3. 

Figure [T4| shows the spectrum of SDSS 
J094603.94+013923.6, an example of a OQ 208-like 
source. Panel (a) shows the continuum decomposition, 
following the same convention as used in the top panel of 
Figure [T] Panel (b) shows the emission-line spectrum after 
continuum subtraction. The red solid line is our Fe II model. 
The blue dashed line marks the position of the peak of Fe II 
A4924 at zero velocity shift. The shift of the Fe II spectrum 
is obvious. Panel (c) illustrates the detailed fitting of the H/3 
and [O III] emission lines. The profile of H/3nc is fixed to 



that of [O 111]. The two blue dashed lines mark the rest-frame 
wavelength of H/3 and [O ill] A5007; both H/3nc and [O ill] 
share the same velocity. A prominent, intermediate-width 
redshifted component is clearly required to fit H/3. Panel 
(d) shows the fit for [O ll] A3727, which also has the same 
redshift as [O ill]. Th is mea ns that this source is not an 
[O III] "blue outUer" (jBgroson 2005). Table |3] fists the 
FWHMs and velocity shifts of Fe 11, [O ill], [O ll], and each 
component of H/3. The velocity shifts are all with respect to 
[O 111]. The width and velocity shift of the redshifted H/3 
component are consistent with those of Fe 11 emission. The 
emission-fine spectrum of SDSS J094603.94H-013923.6 can 
be divided into three systems: (1) H/3, [O III], and [O ll] 
emission lines with narrow (FWHM w 500 km s"') widths 
and no velocity shift; (2) a normal "broad" H/3 component 
with FWHM « 5700 km s"', approximately at rest with 
respect to the narrow lines; and (3) Fe 11 emission and H/3 of 
intermediate width (FWHM « 1500 km s"') redshifted by 
- 1500- 1700 km s-i. 
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TABLE 3 

Emission-line Properties OF SDSS J094603. 94+01 3923.6 





Fe 




H/3 


[O III]" 


[O II] 






Kedshifted 


Broad Nan'ow 






FWHM(kms-') 


1543(55) 


1428(18) 


5730(58) ■■■ I' 


505(6) 


568(18) 


Velocity shift (km s"' ) 


1533(24) 


1756(7) 


373(25) 49(10) 


C 


-22(8) 



Note. — The number in parenthesis is the error. 
" Refers to the line core component used to derive the redshift. 
The width of the H/3nc is fixed to the width of [O III]. 
The redshift of [O III] is used as the systemic redshift. 
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Fig. 13. — H/3 and [O III] region of the composite spectra. The inserted 
plot shows the difference spectra. The vertical dashed line is drawn at 4861 
A. Note the excess on the red wing of H/3. 

OQ 208-like sources can offer an unique view to understand 
the origin of Fe 11 emission and the structure of the BLR. A de- 
tailed study of this class of objects is beyond the scope of this 
paper. We will analyze a sample of such sources discovered 
in our work in a forthcoming publication. For the purposes of 
the present discussion, we simply note: Fe 11 emission is not 
associated with the conventional broad component of H/3 but 
instead originates from the same region emitting the interme- 
diate component of H/3. This is consistent with our findings 
in ^4. II that the width of Fe 11 is systematically narrower than 
that of H/3bc. 

5.2. Where is the Fe II Emission Region? 

This paper has demonstrated that Fe II emission most likely 
does not originate from the same location that produces the 
broad component of H/3. The evidence comes from the sys- 
tematic redshifts and narrower line widths of Fe II compared 
to H/3bc- As discussed below, the simplest interpretation is 
that Fe II emission originates from an inflow that is located 
at the outer parts of the BLR. The H/3bc emission line itself 
shows no systematic velocity shift, and any velocity shift it- 
self is also small. This suggests that H/3bc emission region 
is well virialized and that the width of H/3bc is dominated by 
gravity, making this line suitable for estimating BH masses. 

Our finding that FWHMpe « 3/4 FWHM(H/3bc) suggests 
that the Fe II emission region is farther from the central BH 
than the H/3bc emission region. If the Fe 11 emission region is 



also virialized, so that R oc v"^, then it is about 2 times farther 
from the center than the H/3bc region. On the other hand, the 
systematic redshift of Fe II indicates that the assumption of 
virialization may be incorrect. 

The redward shift of Fe II and the inverse correlation be- 
tween vpe and Lboi/i'Edd favor a scenario in which Fe 11 emis- 
sion emerges from an inflow. To explain the systematic red- 
shift, the inflow on the back side of the accretion disk must 
be obscured by the accretion disk and the torus, so that we 
can only observe the redshifted part on the nearer side. To 
explain the inverse correlation between vpe and Lboi/i'Edd, we 
speculate that the inflow is driven by gravity toward the cen- 
ter and decelerated by the radiation pressure. An increase in 
Lboi/i'Edd would enhance the radiation pressure and lead to a 
decrease of the inward velocity of the inflow. 

Previous studies of the C IV line have shown that it 
tends to be s y stematically blueshift e d ("e.g.. iGaskell 11982 ; 
WilkesI 11984 iMarziani et"an 119961: ISulentic et alJ l2000a : 



Richards et al.l l20"02bh . This finding has led to suggestions 



that high-ionization lines and low-ionization lines originate 
fro m distinct regions (see also the result of reverberation map- 
ping |^terson^&^^ndei[i999|). High-ionization lines such as 
C IV may be emitted from som e kind of outflowing disk wind 
(e.g. Eeighly & Moore 2004; Leighlvll2004 iBaskin & Lao^ 
l2005t fShang et al.l i2007L and references therein), whereas 
low-ionization lines such as H/3 are anchored to a more disk- 
like configuration. This paper adds an additional element 
to this picture. We suggest that in addition to a disk and a 
wind, the BLR has yet another component, one associated 
with inflowing material that produces the Fe II emission. In- 
terestingly, Welsh et al. (2007) recently reanalyzed the spec- 
tral variability of the well-studied Seyfert 1 galaxy NGC 5548 
and suggested, based on the differential lag between the red 
and blue wings of the H/3 profile, that the BLR in this object 
contains an inflowing component. Such an inflow is likely to 
develop from the inner edge of the dusty torus, which may 
connect with the accretion disk. If this picture is correct, the 
redshifted Fe II emission can be used as a probe of the transi- 
tion region from the dusty torus to the BLR or accretion disk. 
More detailed theoretical modelling and observations of the 
Fe II emission region are clearly required to test this picture. 

6. SUMMARY 

Using a large sample of quasars selected from SDSS, we 
have studied the properties of their optical Fe 11 emission, 
especially their velocity profile and velocity shift, with the 
goal of understanding the origin of the nature of the Fe II- 
emitting region. This was accomplished using an improved 
iron template-fitting method, whose reliability has been tested 
using extensive simulations. 
Our findings can be summarized as follows: 
1 . The majority of quasars show redshifted Fe II emission 
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Fig. 14. — Spectrum of SDSS J094603. 94+013923. 6, an example of a OQ 208-like source, (a) Continuum decomposition, (b) Fe model and the observed 
spectrum after subtracting the power law and Balmer continuum. The spectra and the marks plotted in panel (a) and {b) are as the same as those in the top and 
bottom panels of Fig. [T] (c) Fitting of H/3 and [O III] emission lines. The multiple Gaussian components are in blue, and the sum of them is in red. The two blue 
dashed lines mark the rest-frame wavelength of H/3 and [O III] A5007. (d) Fitting of [O II] emission line. The blue dashed line marks the rest-frame wavelength 
of [O II] A3727. Note that the redshifts and widths of [O II], [O III], and H/3nc are almost the same. Note also that we added an additional redshifted Gaussian 
line to fit H/3. The width and velocity shift of this redshifted H/3 component are consistent with those of Fe II emission. See text and Table[3]for more details. 
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with respect to the systemic velocity of narrow-line region (as 
traced by [O III] A5007) or of the conventional broad-line 
region (as traced by H/3). The shift is typically ^ 400 km s"' 
and can be as large as ^ 2000 km s"' . By contrast, neither the 
broad H/3 nor Mg II lines show a systematic velocity shift. 

2. The velocity width of Fe II is systematically nar- 
rower than that of the broad component of H/3. On average, 
FWHMpe « 3/4 FWHM(H/3bc). 

3. The velocity shift of Fe II increases with decreasing Ed- 
dington ratio, decreasing Fe strength, and increasing H/3 line 
width. No clear trends with radio or X-ray properties can be 
discerned. 

4. Composite spectra reveal that objects with large Fe II 
velocity shifts have a tendency to exhibit asymmetric H/3 pro- 
files with an excess red wing. This phenomenon is particu- 
larly notable in a subclass of objects that resembles the pro- 
totype OQ 208. The H/3 profile of OQ 208-like sources show 
a redshifted component of intermediate width that closely re- 
sembles the Fe II emission. 

5. Our results strongly indicate that Fe II emission does 
not originate from the same region of the BLR that produces 
the "traditional" broad component of H/3. Instead, we sug- 
gest that Fe II is associated with the intermediate-width H(3 
component, both tracing an inflowing component of the BLR. 
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